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Introduction 


In  this  paper  a general  formulation  Is  presented  for  the  dynamic 
response  of  a viscoelastic  structural  member  using  a modal  approach. 

This  formulation  will  depend  only  on  the  equations  of  motion  of  the 
structure.  The  formulation  will  apply  to  both  self-adjoint  and 
non-self-ad joint  systems  of  equations  of  motion,  non-homogeneous 
boundary  and  In-span  conditions,  and  arbitrary  excitation  forces  and 
displacements.  In  addition,  proportional  viscous  damping  may  be 
Included  In  a formulation. 

The  first  to  propose  a modal  solution  of  a vlscoelasto-klnetlc 
problem  was  Valanls  [1].  By  assuming  that  Poisson's  ratio  remains 
constant,  the  d3mamlc  problem  can  be  resolved  into  a quasl-statlc 
viscoelastic  problem  and  a dynamic  elastic  problem.  The  latter  problem 
can  be  solved  using  a classical  modal  analysis,  while  the  quasl-statlc 
portion  may  be  solved  using  a correspondence  principle.  Robertson  and 
Thomas  [2]  used  a similar  philosophy,  and  arrived  at  slightly  different 
results  than  Valanls.  These  were  extensions  of  Valanls'  results  Including 
more  general  boundary  conditions  and  non-zero  Initial  conditions. 

Robertson  later  applied  these  results  to  several  specific  members, 
such  as  a circular  viscoelastic  plate  [3],  and  a viscoelastic  beam 
[4].  However,  the  approach  taken  was  not  completely  general,  nor  did 
It  consider  non-self -ad joint  structural  members. 

In  this  paper,  a general  formulation  Is  derived  that  goes  beyond 
the  work  of  the  above  Investigators.  A general  set  of  formulas  Is 
derived  that  explicitly  provides  the  relationship  needed  to  determine  the 
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dynamic  response  of  a viscoelastic  member  due  to  arbitrary  loadings. 

A very  general  set  of  governing  equations  of  motion  can  be  uncoupled 
with  this  approach.  These  equations  may  be  non-self-adjoint  and 
possess  general  non-homogeneous  boundary  and  In-span  conditions.  Two 
illustrations  of  a simple  Volgt-Kelvin  beam  and  an  axisymmetric  viscoelastic 
circular  plate  are  given  to  demonstrate  the  use  of  the  general 
formulation. 
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General  Formulation 


The  dynamic  response  of  a viscoelastic  structural  member  will  be 
formulated  by  using  a modal  analysis.  The  formulation  will  Incorporate 
a damped  viscous  foundation,  which  must  be  proportional  to  the  mass 
and/or  stiffness  of  the  member.  The  differential  equations  governing 
the  linear  motion  of  an  Isotropic  and  homogeneous  viscoelastic  member 
may  be  written  as 


(Xo+Uo)G\^kl  + ^0<^*''l.kk 


32wi  3w^ 
+ c 


with  the  boundary  conditions. 


w^  ■ 0 

and  the  Initial  conditions. 


<#w 

w^(x,0)  - Wj^q(x);  (x,0)  - 


where, 

w^  - displacement  In  the  l*”^  coordinate  direction 

^O’^O  ” constants 

n^  - unit  outward  normal  to  the  surface 

■ prescribed  surface  tractions  on 

”l0'^10  " prescribed  Initial  displacement  and  velocity, 
respectively. 

G ■ time-dependent  relaxation  modulus 
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In  addition,  the  * notation  denotes  the  convolution  Integral, 


f (t)*g(t) 


dx 


(2) 


and  denotes  the  partial  derviative  of  with  respect  to  Note 

that  in  deriving  these  equations,  it  has  been  assumed  that  Polsslon's 
ratio  remains  constant.  As  Valanis  points  out,  several  researchers 
have  shown  that  if  a variation  exists,  it  is  not  substantial,  and  that 
for  incompressible  materials  (or  for  problems  in  which  the  dilation 
is  not  the  primary  mechanism  of  deformation)  this  simplification 
results  in  nearly  exact  solutions. 

The  equation  of  motion  given  by  equations  (1)  may  be  written 
in  the  following  alternate  form: 

* ^ 

G Du  (x,t)  - 2,  A.  (x)3,u(x,t)-F(x,t)  (3a) 

j-1  ^ ^ 

with  the  initial  conditions. 


u(x,0)  - 3u(x,0)  - Uq(x) 

and  the  time-dependent  boundary  and  in-span  conditions. 


(3b) 


G*L^(u,t)  - pj^(x,t)  on  S (3c) 

where  u(x,t)  is  a column  vector  of  dependent  state  variables  including 
displacements  as  well  as  internal  forces.  The  notation  in  these 
equations  is  defined  as  follows: 

F(x,t)  ■ N-dlmenslonal  column  vector  of  body  forces 
D(x),L^(x)  ■ NxN  spatial  matrix  linear  differential  operators 


1.5 


Aj (x)  ■ NxN  spatial  matrix 

P^(x,t)  = N-dimensional  column  vector  containing  the  prescribed 

non-homogeneous  surface  tractions  and  homogeneous  prescribed 
displacements 


Uq(x),Uq(x)  * N-dimensional  vectors  of  prescribed  state  variable 
Initial  conditions 

Note  that  external  proportional  viscous  damping,  such  as  found  In  a 
foundation.  Is  Included  In  the  matrix  A^(x).  In  this  form,  N equations 
of  motion  are  represented  by  a matrix  equation  with  differential 
operator  elements.  It  must  be  kept  in  mind  that  these  equations  only 
apply  to  linear  Isotropic  viscoelastic  materials,  deforming  under 
Isothermal  conditions  with  a constant  Poisson's  ratio.  Nevertheless, 
this  general  form  of  the  equations  of  motion  describes  the  response  of 
a wide  class  of  viscoelastic  structural  members. 

The  dynamic  response  of  a general  viscoelastic  member  may  be  expressed 
In  the  usual  modal  solution  form. 


u(x,t)  - ^ h (t)i()  (x)  (4) 

m-l  ” ” 

In  this  solution,  is  s vector  of  undamped  mode  shapes  corresponding 

to  the  vector  of  dependent  state  variables  u(x,t).  The  coefficient 
h (t)  is  a scalar  factor  which  Indicates  the  relative  contribution  of 

HI 

the  m^^  mode  to  the  overall  response.  In  order  to  find  the  undamped 
mode  shapes,  the  classical  eigenvalue  problem  must  be  solved.  Using 
the  notation  of  equation  (3),  the  undamped  free  motion  problem  takes 
the  form 
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(5a) 


with  the  combined  boundary  and  In-span  conditions. 


L,\p  (x)  * 0 on  S (5b) 

X n 

The  viscoelastic  member  may  be  non-self-adjolnt,  so  that  a blorthogo- 
nallty  relation  will  be  needed  In  deriving  the  elastic  modal  expansion. 
Hence,  use  will  be  made  of  the  adjoint  undamped  free  motion  problem, 

Dii/^(x)  - A2(x)il;^(x)  (6a) 

with  the  boundary  and  In-span  conditions. 


t.'ll  (x)  *0  on  S (6b) 

1 n 

where  the  tildes  refer  to  the  algebraic  adjolnts  of  the  corresponding 
quantities.  The  eigenfunctions,  shown  to 

satisfy  the  blorthogonallty  relation  [5]  given  by 

<>i;„.A2V  - (7a) 

and 

-<(ii^,A2'J^^>  (7b) 

where  6 Is  the  kronecker  delta  and  where  the  notation  <<|)  ,A2<I'  > 
mn  n n 

denotes  the  Inner  product, 

* (A2'l»^)dx  (8) 

Domain 

Solution  of  equations  (5)  and  (6)  can  be  readily  found  for  many 
structural  members.  Thus,  only  the  temporal  coefficients  h (t)  must  be 

ID 
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determined  to  complete  the  modal  solution  In  equation  (4) . To  derive 
an  uncoupled  equation  for  these  coefficients,  begin  with  the  extended 
Green's  Identity, 

<i  , D(G*u)>  - <(G*u),D^  > - B(iii  , (G*u) ) (9) 

“ mm 

The  term  (G*u)  Is  a function  of  both  the  spatial  and  temporal  variables, 

which  Is  assumed  to  be  differentiable  to  the  extent  demanded  by  the 

operator  D.  The  bilinear  form  B(<|»  ,(G*u))  represents  the  non-homogeneous 

m 

boundary  and  In-span  conditions.  It  Is  Identified  by  forming  the  Inner 
product  D(G*u)>  and  then  Integrating  by  parts  with  respect  to  the 

spatial  variables.  Note  that  the  appearance  of  G does  not  effect  the 
Integration  by  parts  because  G Is  a function  of  time.  The  boundary 
and  In-span  conditions  are  then  grouped  to  form  the  bilinear  functional 
B(<(»^,  (G*u))  • Since  the  prescribed  displacements  on  the  boundary  are 
zero,  this  bilinear  form  contains  only  the  nonhomogeneous  surface  tractions 
which  Involve  the  relaxation  modulus  G(t).  In  fact.  It  Is  this  form 
of  (G*u))  that  prohibits  non-homogeneous  displacements  from  being 

prescribed  on  the  boundaries  because  these  may  not  be  expressed  In 
the  form  G*u. 

Before  returning  to  equation  (9) , note  that  a property  of  the 
convolution  Is 

G*(Du(x,t))  - (Dufx,t))*G  (10) 

where  Du(x,t)  may  be  treated  as  an  arbitrary  function  of  x ant  t. 

Since  D Is  a spatial  differential  operator  It  does  not  effect  G(a  function 
of  time)  so  that 
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G*Du  - Du*G  * D(u*G)  - D(G*u)  (11) 

With  this  equation  the  governing  equation  of  motion  (equation  (3a)) 
becomes 

2 

D(G*u)  - I A (x)8  u(x,t)  - F(x,t)  (12) 

j=l  ^ J 

Now  substitute  equations  (6a)  and  (12)  into  the  Green's  identity  in 
equation  (9)  to  obtain 

3<tiJ  ,Aiu>  + .A2U>  + <(G*u),  Aa'i'  > (13) 

m in  m m 

- F>  + 8(1?  , (G*u)) 
m m 

Noting  the  homogeneous  form  of  the  Green's  identity, 

<(G*u),  \z'i>  > * < if  ,A2(G*u)>  (14) 

m m 

and  the  G is  a function  of  time,  it  follows  that 

<if^,A2(G*u)>  “ G*<if^,A2U>  (15) 

or, 

<(G*u)  ,A2'fjjj>  ■ G*<if^,A2U>  (16) 

Using  equation  (16)  in  equation  (13)  gives 

32<if  fA2U>  + 3<if  ,Aiu>  + (X  )^  G*<if  ,A'2U> 
m m m m 

- ,F>  + B(^  ,(G*u)  (17) 

m m 

Invoking  the  proportionality  conditions. 
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(18) 


r 


<\l)  ,A  u>  “ 
m 1 

a common  inner 


(a  + bX  <iii  ,ApU>  ,u) 

lu  in  ^ ni 

product  may  be  factored  out  of  equation  (17)  to  give 


C (t)  + (a  + bX2)c  + (X  )2  G*C 
m m n m in 


(19) 


- <ii)  ,F>  + B(Ti;  ,(G*u))  + b3B(ti;  ,u) 
mm  m 

In  this  equation,  a and  b are  constants  of  proportionality  for  viscous 
damping  such  that 


Aju  = aA2U  - bDu  (20) 

and  C (t)  is  defined  to  be 
m 

C^(t)  = 4^,A2U>  (21) 

so  that  the  generalized  initial  conditions,  5 (0)  and  t (0)  are  found  from  this 

m m 

equation  at  t * 0.  The  coefficients  h^(t)  are  found  by  repeating  the 
steps  outlined  in  Reference  [6]  for  elastic  ntimbers.  Since  the  orthogonality 
relation  remains  the  same  as  for  proportionally  damped  elastic  members, 
it  is  not  surprising  that 


m 


Sn 


(22) 


which  is  identical  to  the  result  previously  obtained. 

Although  a solution  of  the  viscoelastic  dynamic  problem  may  in 
principle  be  given  by  equation  (19),  a more  efficient  form  of  the 
modal  solution  may  be  found  by  extracting  a static  solution  from  the 
modal  expansion.  This  may  be  done  be  rearranging  equation  (19)  into 
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the  form: 


C + (a+bX^)  5 “ H 

m m m m 


where  is  the  generalized  force  term  given  by 


H (t)  =»  JIt  dx  + B(ij>  (G*u))  + b3B(ii;  ,u) 


(23) 


(24) 


Taking  the  convolution  of  both  sides  of  equation  (4)  and  using  equation 
(22)  gives 


G*u  = I 


" (G*C„) 


m 


Q m 


(25) 


m*! 

Now  introduce  equation  (23)  into  (25)  and  after  some  rearrangement  obtain 


oo  l);  H “>  ij; 

ra-1  ''m^  m-1  ^m^ 


C + (a+bX 
m 


m m 


(26) 


The  first  term  in  this  equation  is  the  static  contribution  to  the 
viscoelastic  solution.  As  Valanis  [1]  and  Robertson  and  Thomas  [2] 
have  pointed  out,  a dynamic  viscoelastic  problem  has  a solution  of  the 
form 


u(x,t)  - V^(x,t)  + I 'l'nj(x) 


hfI  ^m 


(27) 


where  v^  is  the  quasi-static  viscoelastic  solution.  Equations  (26) 
and  (27)  are  equivalent  if 
ip  H 


G*v„  - I 


ibfI 


m m 


(28) 
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and 


G*E  - - f ^ r?  + (a+bX2)5  ~| 

m-1  X2  L J 


Then  equation  (26)  becomes 


(29) 


OB  Ip 

G*u  - G*v  + y (G*E  ) -^  (30) 

® n^l  “ ^m 

which  implies  equation  (27).  It  is  now  possible  to  obtain  a more 
thorough  understanding  of  the  terms  in  equation  (30).  The  Green’s 
identity  may  be  expressed  in  the  form 


<il)  ,D(G*v  )>  = <(G*v  ),  DiJ)  > + B(tii  ,(G*v  ))  (31) 

ns  s n ns 

Substitute  equations  (6a)  and  (28)  into  this  equation  to  obtain 

<ip  ,D(G*v  )>  - - I H (t)<t^^,A  * > + B(i  ,G*^^))  (32) 

® m-1  a2q^  “ m 2 n J n s 

where  Integration  of  the  series  is  assumed  to  be  equal  to  a series  of 
the  integrals.  Using  the  biorthogonality  relation  (equation  (7))  and 
in  view  of  equations  (11)  and  (24) , equation  (32)  becomes 


4^,  G*D^g>  - -4^,F>  - b4^.(G*u))  - B(.|»^,(G*Vg))  (33) 

The  bllllnear  forms  in  the  above  equation  contain  contributions  to  the 

generalized  forcing  function  from  the  non-homogeneous  tractions  on  the 

surface.  In  particular,  the  form  G*u))  contains  the  non-homogeneous 

surface  tractions  as  given  by  equation  (3c).  Therefore,  B(\|i  , (G*v  )) 

n s 

implies  that  the  same  boundary  and  in-span  conditions  apply  to  v as 

s 

u when  damping  is  absent  (b~0) . In  addition,  a differential  equation 
valid  throughout  the  volume  is  implied  by  the  remaining  two  term  in 
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equation  (33).  That  Is,  the  following  boundary  value  problem  may  be 
deduced  from  equation  (33) : 

G*D(v^(x,t))  - -F(x,t)  (34a) 

with  boundary  and  In-span  conditions 


G*LiVs(x,t)  - P^(x,t) 


on  S 


(34b) 


The  solution  of  this  boundary  value  problem  (v^(x,t))  will  provide 

the  quasl-statlc  response  of  a viscoelastic  member  Indicated  In  equation 

(27).  Note  that  this  solution  Is  quasl-statlc  In  the  sense  that  v Is 

s 

a function  of  x alone  which  Is  found  for  each  Instant  of  time. 

To  complete  the  d3mamlc  response  solution  In  equation  (27) , an 
uncoupled  equation  for  ^^(t)  must  be  found.  Substitute  this  equation 
Into  the  governing  equation  of  motion  (equation  (3a)  to  give 
/ “ 


G*Dv  + G*D 
s 


I 'J. 

m»l 


Aiv„  + Ai  ( I tj^ 
m“l 


+ AjVg  + A2 


I 


m*l 


mm 


-F 


(35) 


where  It  Is  assumed  that  the  Infinite  series  converges  at  all  points  of 
Interest.  The  first  term  In  equation  (35)  Is  equal  to  -F  which  can 
then  be  cancelled  from  either  side  of  this  equation.  Premultiply  by 
then  Integrate  over  the  volume  to  obtain 

- • " i; 

+ I E^<4;^,Ai^„>  + «iJ„,A2V^> 


m-1 


m 


+ [ <i  ,A2*  > - I (G*E  ) > 

m-1  “ “ “ m-1  “ “ “ 


(36) 
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Note  that  the  validity  of  reversing  the  convolution  and  the  summation 
sign  has  been  assumed.  Apply  equation  (5a)  to  the  above  equation  and 
then  use  the  biorthogonality  relation  for  undamped  members  (equation  (7)) 
to  find 


E (t)  + (a+bx2)S  + x2(G*E  ) - R (t)' 

m m ID  ID  m ID 


(37a) 


where 


(37b) 


and  where  the  proportionality  condition  In  equation  (18)  has  been 

employed.  Note  that  the  condition,  B(<Ji  ,ili  ) 0,  which  implies  that 

n m 

all  boundary  and  In-span  conditions  are  homogeneous  in  the  free  motion 

problem,  has  also  been  used.  The  initial  conditions  for  E (t)  may 

01 

be  found  by  premultiplying  the  Initial  condition  form  of  equation 
(27)  by  '('jjA2  and  then  Integrating  over  x to  give 

« E (0)  _ 

<i  ,A2u(x,0)>  - <!()  ,A2V  (x,0)>  + I (ijj  ,A2<I>  > (38) 

n n s i]i«l  ^ m 

After  applying  blorthogonallty,  the  generalized  initial  conditions 


E^(0)  - j^<*^,A2UQ(x)>  - <<:^,A2V^(x,0)>J 


and  similarly. 


6 (0)  - <tii  ,A2U«(x)>  - <♦  ,A2V  (x,0)>  (40) 

01  01  U 01  8 

An  alternative  form  of  the  right-hand  side  of  equation  (37a)  may  be 
found  which  incorporates  the  complex  compliance  J(t)  of  the  viscoelastic 
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material.  With  this  form,  either  the  experimental  or  modeling  Information 

of  a given  material  may  be  easily  included  in  the  solution  for  E (t) . 

m 

Begin  by  applying  the  proportionality  condition  in  equation  (18)  to 
give 


<ili  .Aiv  > » (a+bX^)  <\l)  .Aov  > - b ,v  ) 
m s mm  s ms 


(41) 


Substitution  of  this  condition  into  equation  (37b)  leads  to 


R - - ^ 
m 


1_ 

Qm 


(c  3+82)  <ii)  ,A2V  > - b B(t(<  ,v  ) 
in  m s ms 


(42) 


where  denotes  the  right-hand  side  of  equation  (37a) , and  c^  is  the 

proportionality  constant  (a+X^b) . The  inner  product  may  now  be  trans- 

m 

formed  by  use  of  the  Green's  identity  into 


<tp  ,A2V  > - <v  ,A2<I>> 
m s s m 


(43) 


From  equation  (6a)  , A2'('  may  be  found  to  be 

m 


^2^ 


Dip 


m 


m 


m 


(44) 


which  when  substituted  into  equation  (43)  gives 


<i|i  ,A2V  > 
m s 


A_  s m 
m 


(45) 


Now  apply  the  extended  Green's  identity. 


8 m ms  ms 


(46) 


to  obtain  from  equation  (45)  the  result. 


m 


m 


(47) 
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The  term  Dv^  may  now  be  transformed  by  taking  the  Laplace  transform  of 
equation  (34a) , 

sGDv  ■ -F  (48) 

s 

where  the  bar  represents  the  corresponding  transformed  variable  and 
G*f  - sGf  (49) 

m 

Equation  (48)  may  be  easily  solved  for  Dv^  to  obtain 

Dv  - - ^ - sJF  (50) 

® sG 

where  the  relation  between  the  Laplace  transfoinns  of  the  complex 
compliance  J and  the  relaxation  modulus  G has  been  used,  l.e. 

— - s j (51) 

sG 

Taking  the  Inverse  Laplace  transform  of  equation  (50)  gives 

DVg  - -J  * F (52) 

Finally,  substitute  this  equation  back  Into  equation  (47),  which  when 
placed  In  equation  (42)  leads  to 


This  Is  the  alternative  form  for  the  right-hand  side  of  equation  (37a) 

which  was  sought.  This  form  of  R (t)  Involves  directly  the  complex 

n 

compliance  J of  the  viscoelastic  material.  When  substituted  back  Into 

equation  (37a)  the  solution  for  E (t)  may  now  be  found  by  taking 

m 


several  approaches.  If  a known  viscoelastic  model  composed  of  springs 


and  dashpots  has  been  selected,  the  compliance  will  be  known  once  the 

model  parameters  have  been  defined.  These  parameters  may  be  prescribed 

a priori  or  chosen  by  fitting  experimental  data  to  the  model  [3],  In 

addition,  experimental  data  may  be  used  directly  In  equation  (53)  without 

any  reliance  on  a particular  model  [4],  In  which  case  equation  (37a) 

would  be  solved  numerically.  When  a model  Is  used,  equation  (37a) 

may  be  solved  by  Laplace  transforms,  which  may  be  Inverted  analytically 

for  simple  models.  An  example  of  how  this  procedure  may  be  carried 

out  for  a Volgt-Kelvln  material  will  be  shown  In  the  following 

section.  Otherwise,  a numerical  Inversion  process  may  be  needed. 

Another  approach  which  may  be  taken  In  finding  the  coordinates  E (t) 

n 

from  equation  (37a)  has  been  suggested  by  Valanls  [1].  When  a Laplace 
transform  of  equation  (37a)  has  been  taken,  a Voltera  Integral  equation 
of  the  second  kind  results.  Then,  any  number  of  techniques  may  be 
used  to  solve  this  Integral  equation. 

To  summarize  briefly,  the  dynamic  response  of  a viscoelastic 
structural  member  on  a proportional  viscous  foundation  may  be  formally 
expressed  by 


» E (t) 

u(x,t)  - v^(*.t)  + I -5 

tBPl  'V 


(54) 


where  'l'jj|(x)  and  are  the  corresponding  undamped  elastic  member 
eigenfunctions  and  norm,  respectively.  The  temporal  coefficients  may 
be  found  by  solving 


E (t)  + (a+bX2)E  + x2(G*E  ) - R (t) 
o HI  ID  m in  n 


(55) 


I 


and  R (t)  may  take  on  either  of  the  two  forms  as  given  by  equations 
m 

(37b)  or  (53).  As  a special  cases,  when  a > b > 0,  a viscoelastic 
member  without  a damped  foundation  may  be  examined,  and  when  the 
bilinear  form  Is  absent  the  member  possesses  homogeneous  boundary  and 
In-span  conditions. 

Illustrative  Examples 

The  formulas  derived  In  the  previous  section  will  be  applied  to 
a Voigt-Kelvln  beam  and  an  axlsymmetrlc  circular  plate.  The  purpose 
In  treating  these  examples  Is  to  demonstrate  the  details  of  the  general 
formulation  for  a simple  member  and  then  to  show  that  the  response  of 
a more  complicated  viscoelastic  member  may  be  uncoupled  by  applying 
the  same  steps  used  for  simpler  members.  In  addition,  the  theoretical 
results  will  be  shown  to  compare  with  those  obtained  by  Robertson  [3], 
whose  less  general  treatment  was  only  valid  for  self-adjoint  circular 
plates  with  no  damping. 

Dynamic  Response  of  a Voigt-Kelvln  Beam 

Consider  the  undamped  Euler-Bemoulli  beam  shown  in  Figure  1. 

It  Is  assumed  that  the  material  may  be  modeled  as  a Voigt-Kelvln 
viscoelastic  material.  The  linearly  varying  distributed  load  q(x,t) 

Is  applied  at  t ■ 0 and  then  maintained  at  that  value  for  all  time. 

That  Is, 

q(x,t)  - j H(t)  (56) 

where  H(t)  Is  the  Heaviside  unit  function,  which  possesses  the  property 

H(t)  * J(t)  - J(t)  (57) 


i 
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The  dynamic  deflection  response  for  this  beam  Is  sought  for  this  loading 
condition. 


The  governing  equations  of  motion  for  a general  viscoelastic  beam 
have  been  given  by  Robertson  [4].  The  equations  may  be  reduced  to 
one  fourth  order  differential  equation  when  shear  deformation  and 
rotary  Inertia  are  neglected.  When  written  In  the  form  of  equations 


the  following  values  are  assigned: 

D - 2(1  * V)  I Aj  - 0 

(58a) 

- -P  F » -q 

(58b) 

u - w(x,t)  Uq  - dp  - 0 

(53c) 

32 

^H^x-0,1  “ ^ ^^2^x-0,l  Tx"^ 

(58d) 

^^l^x-0,1  " ^^2^x-0,£  “ ° 

(58e) 

The  dynamic  response  of  this  viscoelastic  beam  may  be  found  from 
equation  (27)  to  be 


- E (t) 

w(x,t)  - w (x.t)  + I i|»„(x)  (59) 

® mil  Sn  “ 


where  w^(x,t)  Is  the  quasl-statlc  deflection.  For  this  simple  member, 
the  deflected  mode  shapes  f (x)  and  corresponding  eigenvalues  X may  be 


analytically  determined  by  solving  eqiiatlons  (5).  That  Is, 


(x)  - sin  ~ m-  0,+l,+2,... 

HI  % — 


and 


X 

m 


(~)  = 


1.19 


(60) 


(61) 


so  that  the  corresponding  norm  may  be  explicitly  determined  from  equation 
(7b)  to  be 


0^  - Ql/l  (62) 

Note  that  since  the  elastic  beam  is  self-adjoint,  only  the  modes 
are  needed  to  determine  the  response.  Therefore,  only  the  temporal 
coefficient  remains  to  be  determined  in  the  series  portion  of  equation 


(59).  To  find  E (t),  recall  from  equations  (37a)  and  (53)  that 
m 


E (t)  + x2g*E 
m mm 


- — [ / ilJ  (J*q)  dx  ] 

2n  a«-2  n ® 


(63) 


3t 

mm 


where  a * b = 0 since  the  beam  is  undamped.  In  addition,  the  initial 

values,  E (0)  and  fe  (0)  given  by  equations  (39)  and  (40)  are  zero 
m m 

because  the  initial  conditions  Uq  and  Uq  are  zero.  Introduce  equations 
(56),  (60),  (61),  and  (62)  into  the  above  expression  to  arrive  at 


E + X-^G*E 
a in  in 


=iSzl2^]  i(t) 


(64) 


\(mir)5El 

where  equation  (57)  has  been  utilized.  Eliminate  the  time  variable  in 
this  expression  by  taking  the  Laplace  transform  to  obtain 

a 1 + 


s^j 


(65) 


\(iinr)5EI 

where  the  initial  conditions  of  the  complex  compliance  vanish  since  they 

are  proportional  to  E^(0)  2Uid  E^(0)  [37].  Solve  for  the  Laplace  transform 

of  E (t)  (i.e.  E (s))  to  obtain 
n Q 


-•  3 JJL, 


m 


E 

“ (X2  +83j) 


(66) 
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where  for  convenience 


2(-l)V 


“ (mir)  5eI 

The  Laplace  transform  of  the  complex  compliance  for  a Volgt-Kelvln 
material  [7]  Is  given  by 


(67) 


•^0  " 


(k  + s) 


■) 


(68) 


where 


(69) 


and  qQ  and  are  the  spring  and  dashpot  constants  In  the  Volgt-Kelvln 
model.  Introduce  equation  (68)  Into  equation  (66),  and  after  some 
rearrangement  arrive  at  the  expression 


m “ "^m^Ts^^r^TZT^^r^  k ^ (s-r^^)  (s-r2)  ^ 


(70) 


where 


^m*^!  ^m  / 5 , 9 

~2-tr 


(71) 


and  where  the  method  of  partial  fractions  has  been  utilized.  From  a 

table  of  Laplace  transforms,  E (t)  may  be  found  by  Inverting  equation 

in 

(70).  This  results  In 


E^(t)  - -L_[ 

n 


m^  / ' 1 

X _/x^q?-Aq 

m y m l ^o 


, r.t  r-t. 
(r,e  1 - r^e  2 ) 


m 


r-t  r.t  , 

/ 2 1.1  -Kt, 

(e  - e ) - 72  e ] 


(72) 
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From  this  equation  E (t)  may  be  determined.  Returning  to  equation 
(59) , only  the  quasl-statlc  part  of  the  response  remains  to  be  found. 

One  technique  that  may  be  used  Is  the  correspondence  principle  [7]. 

Recall  from  elementary  strength  of  materials  that  the  static  deflection 
of  an  elastic  simply  supported  beam  Is 

w(x)  - [Sx**  - 5x2)i2  + 21**]  (73) 

Using  the  correspondence  principle,  replace  q and  E by  (x/l)  — and  Q(s)/P(s), 

s 

respectively,  to  obtain 

“ ISCTlIq^  ^s(k:  + s)^ 

which  Is  the  Laplace  transform  of  the  quasl-statlc  deflection.  Note 
that  for  a Volgt-Kelvln  solid  [7], 


« q.  (ic  +s)  (75) 

P(s)  ^ 

This  equation  may  be  easily  Inverted  to  give  for  the  quasl-statlc 
deflection  response. 


Wg (x,t) 


. 3x^  - 5x^12  2x1'* 

'■  180  IlqQ 


(1  - 


-Kt. 

e ) 


(76) 


which  decays  as  time  grows  since  k>0.  Hence,  the  solution  of  a Volgt- 
Kelvln  beam  given  by  equation  (59)  has  been  determined.  It  Is  apparent 
that  many  other  methods  of  solution  of  equation  (63)  could  have  been 
pursued.  For  example,  had  J(t)  been  experimentally  determined,  a 
numerical  solution  of  equation  (63)  would  have  been  attempted.  Or, 
this  experimental  data  could  be  fitted  to  a more  complicated  but  more 


j 
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realistic  model.  However,  for  more  complicated  viscoelastic  models  a 


) 


drawback  to  the  previous  method  of  solution  Is  that  an  Inverse  Laplace 
transform  of  a complicated  function  has  to  be  taken.  In  these  cases, 
a numerical  Inverse  Laplace  tr2insform  would  probably  be  warranted. 

Dynamic  Response  of  a Viscoelastic  Circular  Plate 

The  previous  example  was  given  to  demonstrate  the  usage  of  the  general 
formulation  to  find  the  response  of  a simple  member.  The  considerable 
potential  of  the  general  results  will  now  be  exhibited  by  considering 
the  dynamic  axlsymmetrlc  response  of  a viscoelastic  annular  circular 
plate  with  non-homogeneous  moments  and  shears  applied  on  the  boundaries. 

The  equations  of  motion  for  a viscoelastic  plate  have  been  derived 
by  Robertson  13].  In  the  format  of  equation  (3),  the  following  values 
are  assigned: 


^ (i!_  +i  .i«i 

3r2  R 3R  r2  2a 
o 


id  (i_  + 1, 

2R  ''3R  R-' 

o 


( + “ — ) 

^ - R 3R'' 


A,  - phR 

L O 


0 

1 


F - 


u(R,t)  - 


4'(r,t)' 


w(R,t), 


UqCR) 


hK^  i_ 

2R  3R 
o 


1 3 


(77a) 


(77b) 


%«) 


■{o<w' 


(77c) 


The  notation  used  In  these  equations  Is  defined  as  follows: 

<p,w  • slope  and  radial  non-dimensional  deflection,  respectively 

p,  K,  h,R  * mass  density,  shear  area  correction  factor,  thickness, 
and  outer  radius  of  plate,  respectively 
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R * non-dlmenslonallzed  variable  radius,  r/R^  where  r Is 
the  plate  radius 

q » loading  Intensity  on  the  radial  face 

Dq  « hV  12(1  - Vg) 

- h^/  12R  ^ 
o 

The  non-homo geneous  boundary  conditions  are  given  by  equation  (3c). 

That  Is,  on  R “ R^  (the  Inner  radius)  takes  the  form 

(78a) 


and  on  R 


(the  outer  radius) 


(78b) 


Nothing  that  the  bending  moment  along  a circumference  (M)  and  the 
transverse  shear  force  on  a radial  face  (Q)  are  given  by 


" ■ '=*‘’0  < M '■o  R > 

Q . 0,  (,  + ||  ) (80) 


equations  (78)  are  seen  to  represent  an  applied  moment  and  shear  on 
the  Inner  and  outer  plate  radii,  respectively.  These  applied  loads 
are  prescribed  to  be 


using  the  notation  of  equation  (3c) . 
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As  in  the  previous  example,  the  general  formulas  of  the  previous 
section  may  be  used  to  construct  the  modal  solution.  For  this  member, 
the  response  is  given  by  equation  (27)  as 
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Summarv 


In  this  paper  the  dynamic  response  of  a viscoelastic  member  was 
theoretically  formulated.  It  has  been  shown  that  the  uncoupled 
equation  which  results  in  this  modal  approach  is  an  integro-dlf ferential 
equation.  Solution  of  this  equation  is  difficult;  however,  numerical 
solutions  are  possible  [3],  [4],  Note  that  this  paper  advances  the 
past  efforts  of  other  authors  in  that  non-self-ad joint  systems  of  equation 
and  proportional  viscous  damping  may  be  readily  included  in  the  solution. 
Hence,  the  dynamic  response  of  a very  general  viscoelastic  structural 
member  may  be  determined. 
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Figure  1 Voigt-Kelvin  Beam  Under  Variable  Loading 
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NOTATION 


a 


b 

B 


m 


D 

D 

D, 


E 

m 

F 

G(t) 

h 


h 


m 


H(t) 


H 


m 


I 

J(t) 

I 


M 


P p* 

i’  1 


q 


constant  of  proportionality  for  viscous  damping 
N-square  spatial  matrix 

constant  of  proportionality  for  viscous  damping 
bilinear  functional 

viscous  proportionality  constant,  c^  * a + 

N-square  matrix  linear  differential  operator 
N-square  algebraic  adjoint  matrix  differential  operator 

hVl2(l  - Vq) 

modulus  of  elasticity 

temporal  coefficient 

column  vector  of  body  forces 

relazation  modulus  for  viscoelastic  material 

thickness  of  plate 

temporal  coefficient 

Heaviside  unit  function 

generalized  forcing  function  for  structural  member  of 
viscoelastic  material 

moment  of  inertia 

complex  compliance  of  a viscoelastic  material 
length  of  beam 

spatial  matrix  differential  operator 
bending  moment 

unit  outward  normal  to  bounding  surface 

column  vectors  of  non-homogeneous  boundary  and  in-span 
conditions 

applied  loading  intensity 
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Q 


transverse  shear  force  on  a radial  face 


r 

R 

R 

m 

s 

S 

t 

u(x,t) 

V 

s 

V 
w 

X 


classical  norm 
radial  coordinate  direction 
non-dlmenslonal  radial  coordinate  direction 
denotes  right-hand  side  of  viscoelastic  equation 
Laplace  transform  variable 

denotes  bounding  surface  and  In-span  condition  locations 
time 

column  vector  of  dependent  state  variables 

column  vector  denoting  the  quasl-statlc  response  of  a 
viscoelastic  member 

shear  force 

transverse  beam  deflection;  radial  circular  plate  deflection 

displacement  In  the  1^^  coordinate  direction 

general  point  In  the  multidimensional  region; 
one-dlmenslonal  coordinate  direction  In  a rectangular 
coordinate  system 


Greek  Symbols 

h^/12a^ 

6 Kronecker  delta 

mn 

K ^ Volgt-Kelvln  model;  shear  area  correction 

factor 

Xq  constant 

denotes  the  undamped  (or  classical)  frequency  or  eigenvalue 

p constant 

Poisson's  ratio 
o 

temporal  coefficient 
in 
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i 


TT 


n 


3.14159... 
mass  density 
normal  stresses 
dummy  time  variable 
shear  stress 

slope  of  elastic  curve  for  a circular  plate 

column  vector  of  undamped  (or  classical)  mode  shapes  or 
eigenfunctions  corresponding  to  the  nth  frequency 
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